Gauge factor of thick film resistors: outcomes of the variable range hopping model 
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Despite a large amount of data and numerous theoretical proposals, the microscopic mechanism of 
transport in thick film resistors remains unclear. However, recent low temperature measurements 
point toward a possible variable range hopping mechanism of transport. Here we examine how 
such a mechanism affects the gauge factor of thick film resistors. We find that at sufficiently low 
temperatures T, for which the resistivity follows the Mott's law R{T) ~ exp(ro/r)^'''*, the gauge 
factor GF is proportional to (Tq/T)^'^*. Moreover, the inclusion of Coulomb gap effects leads to 
GF ~ {Tq/tY^^ at lower temperatures. In addition, we study a simple model which generalizes the 
variable range hopping mechanism by taking into account the finite mean inter-grain spacing. Our 
results suggest a possible experimental verification of the validity of the variable range hopping in 
thick film resistors. 
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I. INTRODUCTION 

Thick film resistors (TFRs) are composite materials in 
which metallic grains (Ru02, Bi2Ru207, etc.) are em- 
bedded in an insulating glassy matrix. The character- 
istic transport properties of these materials render the 
TFRs particularly suitable as cryogenic thermometers 
and piezoresistive sensors. Despite this successfull prac- 
tical aspect, on the theoretical side the situation is less 
satisfying. The microscopic mechanism of transport in 
TFRs is in fact far from to be understood and several the- 
oretical models have been proposed, none of them how- 
ever being able to provide a completely satisfactory and 
definite answer. 

Among the proposed transport mechanisms in TFRs, 
the most representative ones are those based on tunneling 
within a network of interconnected metallic clusters sepa- 
rated by thin glassy layers hopping between isolated 
metallic grains |^ and phonon-assisted variable-range- 
hopping mechanism |^ . The first model should be consis- 
tent with a highly dishomogenous structure in which the 
metallic phase is organized mainly in large segments sep- 
arated by thin layers of glass while the last two models are 
more suitable for homogeneous dispersions of separated 
metallic grains. The actual situation seems to lay some- 
where in between these two structures. Recent extended 
x-ray absorption fine structure experiments in fact have 
revealed a bimodal distribution of the metallic particle 
sizes in Ru02-based TFRs [Q. According to these mea- 
surements, the metallic phase is organized in large Ru02 
grains or clusters with linear size ranging from ~ 200 
Ato ~ 6000 Aand much smaller Ru02 grains with sizes 
of order 20-70 A. Where the actual current path takes 
place is still under debate. However, It is plausible that 
transport in TFRs takes place mainly along paths con- 
necting the small metallic grains, the large Ru02 clusters 



being in fact too far away (a distance of order 100 Aas 
inferred from transmission electron microscopy Q) from 
each other for an electron to directly hop from a large 
cluster to another large cluster. 

TFRs have also a characteristic temperature depen- 
dence of the resitance R{T). At around room temper- 
ature, R{T) has a minimum and it slightly increases at 
higher T. At lower temperatures, R{T) increases sig- 
nalling a non-metallic behavior. In this low-temperature 
region, R{T) has been observed to follow a exp(To/T)^ 
behavior with x ~ 1/4 m or a; ~ 1/2 [§ depending 
on the room temperature resistance value of the sam- 
ples. Recent measurements ^ have shown that R{T) ~ 
exp(To/T)^ with x — 1/4 with a cross-over to a; = 1/2 
as the temperature is lowered for the most resistive sam- 
ples. Such kind of temperature dependence is typical of 
Mott variable-range-hopping (VRH) mechanism of trans- 
port Q affected by Coulomb interaction between grains 
(Coulomb gap effect 1^). Hence, the transport proper- 
ties of TFRs at low temperatures seem to be qualitatively 
described by ordinary VRH as other disordered systems 
like n-type CdSe ||l^. However, the values of Tq and Tq 
that best fit the experimental data would lead to values 
of the optimal hopping length Rh of the same order of 
the localization length ^, whereas the VRH mechanism 
should be observed only for i?^ ^ f The validity of 
the VRH mechanism is therefore not clear and additional 
informations are required to analyse the validity of the 
VRH approach in describing transport in TFRs. 

Here we propose that a quantity which can be helpful 
and easily measurable is the strain sensitivity, or gauge 
factor, and its temperature dependence. The gauge fac- 
tor (GF) relates the variation of the total resistance R{T) 
with an applied strain e and it is defined as follows: 



GF 



5R{T) 

Tmr)'' 



(1) 



where 5R{T) is the variation of R{T) under the applied 
strain s. TFRs have high values of GF, typically between 
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2 and 35 at room temperature , and it is precisely this 
property that renders these systems particularly suitable 
for piezoresistive sensor applications. Here we are inter- 
ested in studying how the VRH mechanism affects GF 
and its temperature dependence. We find that if the low 
temperature transport is due to the VRH mechanism then 
the gauge factor acquires a characteristic temperature de- 
pendence which could be experimentally determined. 

In the next section we use the asymptotic regimes pre- 
dicted by the VRH model to extract qualitative behaviors 



for GF. In Sec. Ill , a simple model is introduced with the 



aim of studying the effect of finite inter-grain distances 
and to test the validity of the qualitative behaviors of 
GF. 



II. QUALITATIVE RESULTS 

In this section we analyze the transport properties of 
TFRs from the point of view of VRH model. To this end, 
we must assume drastic simplifications and we model 
therefore the TFRs structure as given by an homoge- 
neous dispersion of small metallic grains (of typical sizes 
20-70 A) in the amorphous matrix and we neglect any 
possible effect of the larger metallic clusters on the bulk 
resistivity. As we discuss below, the assumption of ho- 
mogeneity probably leads to a strong under-estimation 
of the room temperature GF value. However, instead of 
absolute values, here we are mainly interested in relative 
changes of GF as the temperature is varied. 

An additional simplification we assume is to disregard 
the high temperature rise of R(T) beyond the temper- 
ature of minimum. Although the origin of this feature 
is still debated, it should nevertheless been given by 
electron-phonon intra-grain scattering, thermal expan- 
sion effects or maybe a combination of both. Within 
these approximations, at sufficiently high temperatures 
transport is governed by tunneling between adjacent 
grains which are separated by a mean distance d. In 
this high temperature regime therefore 



R{T) = i?o oc exp(2d/0 



(2) 



where ^ is the localization length and it is related to 
the barrier potential V between adjacent grains as ^ = 
h/\/2mV, where m is the electron mass. An applied 
strain modifies d and the resulting longitudinal GF be- 
comes 



GF = GFo = 2d/^. 



(3) 



In the above expression it is implicitly assumed that the 
electrons hop along paths mostly parallel to the direction 
in which strain is applied. However, the actual micro- 
scopic paths are the result of hopping processes also along 
directions perpendicular to the applied strain. In this 
way, the longitudinal gauge factor is something less than 



2d/^ and the amount of reduction depends on specific 
material properties. However, as shown in appendix 
such a reduction can be estimated by introducing a phe- 
nomenological parameter x which measures the percent- 
age of hops along directions perpendicular to the strain, 
leading to the following expression for the longitudinal 
gauge factor: 



GFo = i^Mi^2d/e, 



(4) 



where v is the Poisson ratio which varies between 0.2 
and 0.4 for typical TFRs As shown in appendix 
< X < 1/2 and consequently 0.1 (2(i/^) < GFq < 2d/C 
Since ^ is of order of d, these values of GFq are much less 
than those measured in TFRs (2 < GFq < 35). We think 
that a strong enhancement of GFq could be achieved 
when the assumption of homogeneity is relaxed. In fact, 
a strongly non-homogeneous distribution of metal grains 
in the glassy matrix could lead to a local strain distribu- 
tion very different from the averaged one. In particular, 
the paths along which the current flows can enter regions 
of concentrated strain leading therefore to an enhanced 
total GFq. Studies along this direction are currently un- 
der developement. 

Having obtained the high temperature resistance and 
strain sensitivity, we analize now the situation at low 
temperatures. According to the VRH model, at sufh- 
ciently low temperatures hopping to adjacent grains is 
no longer favourable and R(T) acquires a temperature 
dependence. The question is whether the VRH mecha- 
nism of transport affects the low temperature regime of 
GF as well and if this variation can be experimentally 
tested. 

In the VRH model, the transport properties are gov- 
erned by the probability Pij that an electron hops from 
grain i to grain j [S| : 



cx exp 



-2r.,/e- 



(5) 



where is the distance between grains i and j and 
E{rij) is the energy threshold the electron experiences 
in hopping the distance . Assuming that a finite frac- 
tion of the grains is charged, the energy E{rij) is made 
of two contributions [llOi: 



Eir,,) = ^ + ^, 
grf. nr.j 



(6) 



where g 
charge, 
sionless 
i), the 

a grain a distance r.y . This term is proportional to r-^-^ 
because the probability an electron at site i has to find 
a site j with energy much nearer to its own scales as rf - . 



is the electron density of states, e is the electron 
K is the dielectric constant, f3 and 7 are dimen- 
constants. In the second right hand of equation 
first term describes the energy needed to hop to 
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The second term in Eq. (|^) describes the Coulomb en- 
ergy due to finite charging of the grains, as it is expected 
by the fluctuations in the energy levels in the grains. | |T^ ] 
This term is responsible for the opening of the Coulomb 
gap. @ 

At sufficiently low temperatures, the resistance is gov- 
erned by 1/Pmax where Pmax is the maximum of Pij. 
Hence, the optimization of the exponential in the hop- 
ping probability (^) leads to a temperature dependence 
of R{T) characterized by two cross-over temperatures: 



2048^ 



(7) 



which stems from the first term in the second hand side 
of Eq.d), and 



T' 



876^ 



(8) 



which is given by the Coulomb interaction. Usually Tg < 
To and i?(T) ~ e^p{To/T)^^^ for < T < Tq while 
i?(T) ~ exp(T(;/T)i/2 for T < T^. 

Depending on the particular temperature range, the 
strain sensitivity of R{T) is governed by the strain sensi- 
tivity of To or Tq. Let us first focus on the temperature 
range in which R{T) ~ exp(To/r)^/''. An applied strain 
e leads to following relation between the resistance vari- 
ation 6R and STq: 



SR 
~R 



4 I T 



To 



(9) 



Tq is inversely proportional to the electron density of 
states g which is proportional to the density of grains 
n. Assuming that an applied strain affects the density 
of metallic grains leaving the volume of the grains un- 
changed we obtain STq/Tq — —Sg/g = e{l — 2v) and 
the gauge factor GF = SR/eR becomes: 



GF = 



1 - 2i/ 



1/4 



(10) 



GF therefore increases as T ^Z"* as the temperature is 
lowered provided T is in the range for which R{T) ~ 
exp(To/r)i/'' holds true. 

By further lowering the temperature, transport be- 
comes affected by the Coulomb interaction and R{T) 
crosses towards the ~ exp(TQ/T)^/'^ regime. Hence, un- 
der applied strain, SR/R reduces to: 



6R 



1/2 



6n 



(11) 



The variation 5Tq is driven by the strain dependence of 
the total dielectric constant k and, in full generality, it 
can be expressed as: 



T' 



Sk , , 
= Ae 1 

K 



2iy) 



(12) 



where A is a dimensionless parameter and its value de- 
pends on the dielectric constants of the metallic and 
glassy phases. In appendix |b| we provide an explicit 
expression of A based on the Maxwell-Garnett formula 
for the dielectric response of small metallic particles sus- 
pended in a dielectric. By using (|ll] ) and (|l^), the gauge 
factor in the Coulomb regime reduces to: 



GF 



1/2 



(13) 



The Coulomb effect is therefore reflected in a cross-over 
from a T~^/'* to a T~^/^ dependence of GF as the tem- 
perature is sufficiently lowered. 



III. SIMPLE MODEL FOR THE GF 
TEMPERATURE DEPENDENCE 

Summarizing the main results obtained in the previous 
section, the temperature dependence of the resistance 
and the intrinsic GF as T is lowered is characterized, 
within the VRH theory, by three well distinguishable 
trends: 



R ~ exp(2d/0 ^ GF - 2d/^ 
R ~ exp 

R ~ exp 



,1/4 
1/2 



(14) 



Note that for each behavior listed above, GF is pro- 
portional to Inp(r). The last two behaviors originate 
from the optimization of the inter-grain hopping prob- 
ability Pij of Eq.(||). It is possible in fact to define a 
universal functional form for i?(T) which contains both 
exp(ro/T)i/'' and exp(r^/r)i/2 as limiting values and 
which fits well transport data of Indium doped CdSe sam- 
ples. ||l^ However, for sufficiently high temperatures, the 
optimization of Pij requires optimal hopping distances of 
order or less than the mean distance d between adjacent 
grains. It would therefore be more correct to formulate 
the transport problem in such a way that the value of d 
is explicitly included. Here we provide a simple version 
of such a formulation capable of describing automatically 
the limiting behaviors listed in Eq.(|l^. 

A simple general equation for the conductance G can 
be defined as follows. A single hopping process between 
grains at distance and with a hopping probability Pij 
can be regarded as a resistive element with resistance Rij 
proportional to 1/ Pij . Hence, the total resistance can be 
constructed by interpreting paths characterized by differ- 
ent hopping distances as resistors in a parallel geometry. 
The resulting conductance should therefore be given by 
a weighted summation of the hopping probabilities Pij 
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for all values of r^- > d. In this way, by employing a con- 
tinuous approximation, we define the conductance G{T) 
as: 



G{T) = a / drf{r) exp 



E{r) 



KbT 



(15) 



where E{r) is given in Eq.® and a is a constant intro- 
duced for dimensional purposes. In the above equation, 
/(r) is a weight function which takes into account the 
different path geometries. For our purposes, the detailed 
structure of /(r) is unnecessary and we retain only its 
main feature which is given by a lower cut-off represent- 
ing the mean distance d between adjacent grains. Hence, 
we approximate f{r) by the step function 9{r — d). In 
this way, equation ( p5| ) reduces to: 

f+°° ( E(r) \ 

G{T) =aj^ dr exp (^-2r/e - j . (16) 

For temperatures larger than 



2dKY> 



(17) 



equation (|l^) gives the proper high temperature limit 
(^, while for T ^ Ts the integral can be estimated by 
optimizing the exponential, leading therefore to the VRH 
regime. Concerning the gauge factor resulting from (16), 
it should be noted that an applied strain e affects both 
the interval of integration and the energy threshold func- 
tion E{r) through the electron density of states g and 
the dielectric constant n. Therefore, by using the strain 
derivatives of g and n introduced in the previous section, 
GF can be expressed as the sum of three different con- 
tributions, GF = GFd + GYg -f GF« where: 



GFd = 



2d 



i-x/(i-^) 



1 



G{T) _E 
Go 



{d)/KuT 



(18) 



GF 

^ K^T G{T) 



+ 00 



GF, = 



l-2v \a 
KbT G{T) 



dr-^ exp I —2r/^ 



dr exp I — 2r/^ 

KT 



E{r) 



KbT 



(19) 



E{r) 



KbT 



(20) 



In equation ( |iq ) Go = exp(— 2(i/^)/2 is the asymptotic 
conductance in the high temperatures limit. 

In Fig. la we plot the resistance R{T) obtained nu- 
merically from Eq.(|l^) in the absence of Coulomb re- 
pulsion (7 = 0) for 2d/^ = 0.2, 2, and 5. The resis- 
tance data are in units of the high temperature limit 
i?o — 1/Go and the temperature is given in units of 
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FIG. 1. (a): temperature dependence of the resistance cal- 
culated from Eq. (Ilq) without Coulomb interaction (7 = 0) 
and for the different values of 2d/5 reported in the parenthe- 
ses (solid lines). The dashed lines fit the numerical data with 
exp(fo/r)^/* where Tb ~ 0.72To. (b): corresponding tem- 
perature dependence of the gauge factor for the same values 
of 2d/5 as in Fig. la and for v — 0.3 and x = 1/4. 

To [see Eq.(0)]. The curves show two distinct regimes 
as a function of temperature due to the finite values of 
the averaged inter-grain distance d. By using Eq.(^7|), 
the cross-over between these two regimes takes place 
roughly at T - = 0.105To/(2d/^)4, so that as 2d/^ 
increases the cross-over temprature is reduced. For 
T > Ts [(To/T)i/4 ^ l.755(2d/^)] the resistance shows 
an activated regime while at sufficiently low tempera- 
tures [(To/T)i/4 > l.755(2d/C)] it becomes of VRH type. 
The latter behavior is signalled by the dashed lines which 
fit the resistance with R{T)/Rq cx exp(ro/r)^/'', where 
To ~ O.TSTq- The gauge factors GF associated to the 
resistance curves of Fig. la are plotted in Fig. lb for 
V = 0.3 and x = 1/4. A 2c?/^ dependent cross-over is ob- 
tained for GF as well. For temperatures higher than T,, 
GF is independent of T and approaches the limit GFo 
given in Eq.(^). On the contrary, by lowering T, GF 
acquires a temperature dependence and for T ^ T, in- 
creases as 1/T^/'^ in agreement with Eq.(p^). From the 
results of Fig. lb, the 1/T^/^ dependence of the gauge 
factor holds true in a wide range of temperatures only 
for quite small values of GFo. 

The effect of Coulomb repulsion on R{T) and GF is 
displayed in Fig. 2 where the case 2d/^ = 2 is investi- 
gated for different values of the Coulomb parameter Tq 
defined in Eq.(|). In Fig. 2a, for ^/Tq ^ the resis- 
tance approach the low temperature asymptotic regime 
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FIG. 2. (a): temperature dependence of the resistance with 
the inclusion of Coulomb interaction for 2d/^ — 2 and Tg/To 
as reported in the parentheses (solid lines). The dashed lines 
fit the numerical data with exp(fo/r)^''^ where Tq ~ O.gTo. 
(b): temperature dependence of the gauge factor for the same 
values of Tq/Tq as in Fig. la, i/ — 0.3, x = 1/4, and A = 5 
(solid lines). The dashed lines fit the low temperature regime 
of GF with {f^/Ty^'^. 

R{T)/Ro oc exp(f^/T)i/2 with fj^ ~ Q.<dT^ (dashed 
hnes). The Coulomb effect is reflected also in the temper- 
ature dependence of GF, Fig. 2b, where the curves have 
been calculated by using X — 1, v — 0.3 and x = 1/4. 
The set of curves with Tq 7^ show the low temperature 
asymptotic behavior GF oc (Tq/T)^/^ with values of Tq 
somewhat larger than those obtained from the resistivity 
data (Tq ~ I.OSTq). Finally the temperature dependence 
of GF for different values of A is plotted in Fig. 3 for the 
case 2d/£, = 2 and T(^/To = 0.05. 

IV. DISCUSSION AND CONCLUSIONS 

The numerical solutions of equations (|l6|)-(p0|) confirm 
the qualitative results given in the previous section, i.e., 
as the resistance crosses over the VRH regime, the gauge 
factor acquires a temperature dependence which is qual- 
itatively summar ized in Eq.(|l^. 

The model here presented is very simplified and sus- 
ceptible of various improvements. For example, from 
statistical arguments, it would be certainly more correct 
to evaluate the temperature dependence of K{T) within 
the percolation theory for transport |jl^. This could be 
achieved by explicitly taking into account in the eval- 
uation of the critical path the lower cut-off d for the 
hopping distance and the strain dependence of g and n. 
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FIG. 3. Gauge factor calculated for To/To = 0.05, u = 0.3, 
X = 1/4 and for different values of A (reported in the paren- 
theses) . 

As already stressed before, an additional simplifying as- 
sumption we have used is that of homogeneity. However 
the structure of TFRs is very complex, being made by a 
glass with dispersed metallic grains with bi-modal distri- 
bution of sizes. In this situation, under an applied strain, 
the stress distribution inside the sample could strongly 
modify GF as compared to the homogeneous case and to 
settle quantitatively this modification is certainly a very 
important issue. 

A final open question regards the general applicability 
of the VRH concept when electrons hop between grains 
of finite size. In fact, in reality, the discrete electron en- 
ergy levels of the grains are smeared by finite life-time ef- 
fects given by the electron-phonon and electron-electron 
interactions together with the inter-grain tunneling cou- 
pling. At sufficiently low temperatures, the life-time ef- 
fects are small and the discreteness of the energy levels 
requires phonon assisted hopping for electron transport. 
However, at sufficiently high temperatures, the energy 
smearing due to the finite electron life-time can be of the 
order of the level spacing. In this situation it can be 
easily shown that electrons can directly tunnel between 
grains without the assistance of phonons. Of course, a 
crucial role is played by the grain size which governs the 
energy level spacing. To our knowledge, the problem of 
how the VRH between grains is affected by the inclusion 
of finite life-time effects has never been addressed and its 
investigation is certainly of the most valuable interest. 

In conclusion, the results here reported provide a clear 
prediction of the temperature dependence of the gauge 
factor in the hypothesis that transport is well described 
by the classical VRH theory. Hence, the simultaneous 
detection of the temperature dependences of the sheet re- 
sistance and the gauge factor can provide an useful exper- 
imental tool to test the validity of the VRH mechanism 
in TFRs. We note that some experimental data show 
that GF for commercial TFRs is largely independent of 
T from room down to cryogenics temperatures ||l^ . Our 
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analysis would indicate therefore that these experimen- 
tal results are not compatible with the VRH mechanism 
of transport. However, it is important to stress that the 
VRH model predicts a strong temperature dependence of 
GF only when the resistance crosses over a exp(To/T)^ 
type of regime {x — 1/2 or 1/4), a test that has not been 
performed in Ref. Instead, a firm experimental test 
requires measurements of the temperature dependences 
of both R{T) and GF. 
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APPENDIX A: 

In this appendix we estimate the reduction of the gauge 
factor due to the contributions of hopping processes per- 
pendicular to the direction on which strain is applied. 
To this end, we first consider an operative definition of 
the gauge factor. Let us consider a cantilever beam in 
the x-y plane with axis lying on the x direction on top 
of which is deposited the TFR. The thicknesses of the 
cantilever bean and of the TFR are measured in the z 
direction. According to elasticity theory, bending of the 
cantilever beam produces a finite strain Sxx = e in the 
X direction and no strain in y (syy = 0). Moreover, if 
is the Poisson ratio, the strain of the TFR along the z 
direction is Szz — —^vj^l — v). 

To measure the longitudinal GF of the TFR, a po- 
tential difference must be applied along the x direction. 
If the microscopic paths are made of discrete hopping 
processes exclusively along the direction of the field, the 
longitudinal GF would be equal to 2d/^ just as in Eq.(|[). 
Conversely, the transversal GF would be zero because in 
this case the field is applied in the y direction for which 
Eyy — 0. Let us consider now the situation in which 
there is a finite percentage x of hops along directions 
perpendicular to that of the applied field. The resulting 
resistance Rx , when the field is applied along the x direc- 
tion, becomes Rx = {1 — x)Rx + x{Ry + Rz), where Rx, 
Ry and Rz are resistances for hops in the x, y and z di- 
rections, respectively. Although Rx, Ry and Rz have the 
same average value, they differ when strain is applied. In 
fact: 



5Rx 



2d SRy 



0; 



SRz 



2d 



Hence the longitudinal gauge factor GF^ becomes: 



(Al) 



GFr 



SRx ^ (1 - x)SRx + xSRz 
ZxxRx exx[{i-x)Rx+x{Rv + Rz)] 

\-x/(l~v) 2d 



1 



X 



(A2) 



By following the same reasoning, the transversal gauge 
factor G¥t can be estimated by realizing that the re- 
sistance for a field applied in the y direction is Ry — 
(1 — x)Ry + x{Rx + Rz) and consequentely: 



GFt = 



5Ry 

^xxRy 



x{5Rx + 5Rz) 



exx[{i - x)Ry + X{.Rx + Rz)] 
X[l - - ly)] 2d 



1 + X 



(A3) 



If the Poisson ratio v is known, the phenomenological 
parameter x can be estimated by measuring the ratio 
between the longitudinal and the transversal GFs which, 
according to ([A^) and (AS), is given by: 



1 



x-^ 



GF/ 



GFt x(1 - 2!^) 



(A4) 



Experimentally this ratio is found to be larger than the 
unity and consequently the above expression leads to x < 
1/2. 



APPENDIX B: 

Here, we evaluate the strain dependence of the dielec- 
tric constant for TFRs. We model the structure as given 
by spherical metallic particles dispersed in a glassy ma- 
trix. The metallic and glassy phases occupy volumes 
Vin and Vg, respectively, so that the volume fraction of 
the metallic phase is (/) = T4i/(Kn + V^) ~ Vra/Vg for 
low metallic concentrations. Within the dipole approx- 
imation, the total dielectric constant k is given by the 
Maxwell-Garnett formula [1T5|,E6|: 



1 



(Bl) 



where Kg is the dielectric constant of the glass. For suf- 
ficiently small values of (p, the variation 5k, given by an 
applied strain e is: 



5k, 

5k — K — - — 3Kg(/)e(l — 2v), 

Kg- 



(B2) 



where we have considered the metallic particles as per- 
fectly rigid bodies compared to the glass so that 5(1) — 
-4>e{l-2v) m. 

To find the variation ^Kg we should know an ex- 
plicit expression for Kg which is however a difficult prob- 
lem. Nevertheless, experiments suggest that dielectric re- 
sponse of glasses is under several aspects quite similar to 
that of ionic solids and, in particular, it has been shown 
that the variation of the glass polarizability with pres- 
sure resembles that of the ionic crystalline compounds 
p6[ . By starting with this latter observation we could 
argue that, for our purposes, the dielectric constant for 
simple glasses is roughly given by: 
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1 + 87ra/3fi 
1 -47ra/317' 



(B3) 



where fl has dimension of volume and a ~ 1/K is the 
ionic polarization wit h K being the elastic constant. An 
expression similar to ( B3 ) is also recovered in theoretical 
analyses based on coarse-graining approach (there, 
ri is a course-graining volume). An applied strain modi- 
fies Kg via the volume J7 and the elastic constant K. The 
strain dependence of the latter quantity can be guessed 
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tains K = l^on(''o) ~ l/''o- This result suggests that a/fi 
scales as 1/ [distance] and consequently from Eq.(B3) one 
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The above expression plugged into Eq.(B2) leads finally 
to: 
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where the last equality defines the dimensionless quantity 
A used in the main text. 



[1] G. E. Pike and C. H. Seager, J. Appl. Phys. 48, 5152 
(1977). 

[2] P. Sheng and J. Klafter, Phys. Rev. b 27, 2583 (1983). 

[3] N F. Mott, J. Non-Cryst. Solids, 1, 1 (1968). 

[4] C. Meneghini, S. Mobilio, F. Pivetti, L Selmi, M. Pru- 
denziati, and B. Morten, J. Appl. Phys. 86, 3590 (1999). 

[5] S. Piccinini, Ph. D. thesis. University of Modena, 1998. 

[6] Q. Li, C. H. Watson, R. G. Goodrich, D. G. Haase, and 
H. Lukefahr, Cryogenics 26, 467 (1986); B. Neppert and 
P. Esquinazi, Cryogenics 36, 231 (1996). 

[7] W. Schoebe, Physica B 165-166, 299 (1990). 

[8] M. Affronte, M. Campani, S. Piccinini, M. Tamborin, B. 
Morten, and M. Prudenziati, J. Low Temp. Phys. 109, 
461 (1997). 

[9] A. L. Efros and B. L Shklovskii, J. Phys. C 8, 249 (1975). 

[10] A. Aharony, Y. Zhang, ynd M. P. Sarachik, Phys. Rev. 
Lett. 68, 3900 (1992). 

[11] M. Prudenziati, Handbook of Sensors and Actuators (El- 
sevier, Amsterdam, 1994), p. 189. 

[12] E. Cuevas, M. Ortuno, and J. Ruiz, Phys. Rev. Lett. 71, 
1871 (1993). 

[13] We are not aware of measurements reporting the Young 
modulus E of bulk Ru02. However, for metallic ruthe- 
nium the Young modulus is quite large {E ~ 450 GPa) 



7 



